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I.  FINAL  PROGRESS  REPORT  FOR  GRANT  NUMBER  W911NF-09-1-0347 

The  aim  of  this  project  is  to  explore  the  strongly  correlated  states,  more  specifically  Bose-Einstein  condensation  of 
indirectly  bound  excitons,  in  bilayer  graphene.  During  this  final  reporting  period,  we  have  continued  our  work  on 
the  extention  of  the  many-body  calculations  to  an  atomistic  model,  examined  the  electron-phonon  coupling,  studied 
the  possibility  of  finding  the  same  type  of  condensation  in  a  3D  topological  insulator,  and  started  our  work  on  a 
time-dependent  theory  to  understand  the  behavior  of  the  time-dependent  many-body  systems.  In  this  final  progress 
report,  we  will  outline  our  most  important  results  and  draw  some  final  conclusions  from  our  3  year  project. 


II.  OVERALL  SCIENTIFIC  ACCOMPLISHMENTS 
A.  Overall  Scientific  Accomplishments:  Electron-Phonon  Interactions  in  Dipolar  Superfluids 

In  Fig.  1,  we  show  a  schematic  of  the  graphene  pseudospin  (G-PsS)  device  we  consider  to  understand  the  role 
of  the  electron-phonon  interaction  in  multi-layer  graphene  structures.  The  proposed  device  is  comprised  of  two 
monolayers  of  graphene  with  lengths  of  30  nm  and  widths  of  10  nm.  The  layers  are  considered  to  be  perfectly 
registered  ideal  sheets  of  zigzag  graphene,  such  that  the  lattice  sites  of  each  layer  are  aligned.  Each  monolayer  is 
terminated  with  an  independent  set  of  metal  contacts  from  which  current  may  be  injected  and  extracted.  We  assume 
an  ideal  model  for  these  contacts,  which  do  not  cause  unintentional  doping  of  the  graphene  layers1 .  The  graphene 
monolayers  are  separated  by  a  1  nm  layer  of  Si02  which  serves  as  a  tunnel  dielectric  with  a  corresponding  single 
particle  tunneling  amplitude  of  Asas  =  1  peV.  Each  layer  is  individually  gated  to  independently  control  the  layer 
carrier  concentrations  with  the  gates  separated  from  the  layers  by  20  nm  of  Si02 .  In  order  to  most  effectively  drive 
the  interlayer  condensate  current,  we  apply  the  bias  voltage  in  the  drag-counterflow  configuration2  in  which  positive 
interlayer  voltage  is  applied  to  the  top  left  contact  ( Vtl )  and  an  equal  negative  interlayer  bias  is  applied  to  the  top 
right  contact  (Vtr)  while  the  contacts  in  the  bottom  layer  remain  grounded  (Vbl  =  Vbr  =  0  V).  The  interlayer 
current  discussed  here  is  then  measured  flowing  from  the  top  left  contact  to  the  bottom  left  contact. 

To  model  our  system,  we  use  a  tight-binding  Hamiltonian  to  represent  each  of  the  layers 

Slayer  =  £  T  |  i)(j  \  +Vt  \  i)(i  |,  (1) 

<i,3> 

where  lattice  points  i  and  j  are  first  nearest  neighbors,  r  =  — 3.03eV  is  the  nearest  neighbor  hopping  energy  for 
the  pz  orbital  of  graphene,  which  allows  for  the  unique  low-energy  linear  dispersion  at  the  K  and  K'  points  in  the 


FIG.  1.  Schematic  illustration  of  the  device  under  consideration  in  this  work.  Here  we  have  two  monolayers  of  graphene  separated 
by  1  nm  Si02-  Each  layer  has  a  separate  set  of  contacts  (Vtl  and  Vjr  for  the  top  layer  and  Vbl  and  Vbr  for  the  bottom  layer) 
from  which  current  may  be  injected  and  extracted.  Top  and  bottom  gates,  Vtg  and  Vbg,  are  used  to  adjust  the  layer  carrier 
concentrations. 
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Interlayer  Bias  (V) 


FIG.  2.  Current- Voltage  characteristic  obtained  from  NEGF  simulations  comparing  ballistic  (dotted  lines)  and  phonon  cases  (solid 
lines)  for  two  gate  voltages,  Vb=0.4V  (black)  and  Vg=0.2V  (red). 


Brillouin  zone.  The  full  system  Hamiltonian  is  formed  as 

+  ^2  (i 

/i=x,y,z 


'Hsys  — 


H-layer  0 
0  H layer 


(2) 


with  the  interlayer  interactions  including  both  single  particle  tunneling  and  the  mean-field  many-body  contribution, 
A,  coupling  the  two  layers  using  a  local  density  approximation,  which  drives  the  phase  transition7 .  In  Eq.  (2),  p 
represents  a  vector  that  isolates  each  of  the  Cartesian  components  of  the  pairing  vector,  and  represents  the  Pauli 
spin  matrices  in  each  of  the  three  spatial  directions.  The  on-site  potential  energy  VL  is  calculated  via  a  3-dimensional 
Poisson  solver.  The  transport  quantities  are  calculated  via  the  non-equilibrium  Green's  function  (NEGF)  formalism3 
self-consistently  with  the  interlayer  interactions  and  the  electrostatics.  We  introduce  carrier-phonon  scattering  to 
our  Hamiltonian  within  the  self-consistent  Born  approximation.  We  include  the  relevant  self-energy  terms  for  both 
elastic  acoustic  phonon  scattering  (with  deformation  potential  Dac  =  19  eV)4  and  inelastic  optical  phonon  scattering 
processes  (with  optical  deformation  potential  Dop  =  lxlO9  eV/cm  and  energy  hujo  =  164.6  meV)5,  in  which  both 
interband  and  intraband  scattering  are  included  by  utilizing  a  real-space  basis6.  We  have  set  the  system  temperature 
to  be  Tsys  =  300  K  for  all  calculations. 

In  G-PsS  devices,  the  most  important  parameter  to  device  performance  is  the  critical  current  (7C),  or  the  maximum 
interlayer  current  that  the  system  can  sustain  by  adjusting  the  interlayer  phase  relationship.  Ic  is  defined  as7 


qWps 

hLc 


(3) 


In  Eq.  (3),  q  is  electronic  charge,  W  is  the  width  of  the  system,  ps  is  the  pseudospin  phase  stiffness,  and  Lc  is  the 
coherence  length.  The  coherence  length  in  pseudospin  systems  expected  to  form  a  condensate  phase,  defines  the 
maximum  length  which  quasiparticles  propagate  in  the  system  before  undergoing  retroreflection.  When  a  quasi¬ 
particle  is  retroreflected,  the  system  launches  an  exciton  to  conserve  current7 .  The  process  is  similar  to  Andreev 
reflections  in  a  superconductor,  except  that  here  electrons  are  injected  in  one  layer  and  are  reflected  as  electrons  with 
opposite  momentum  in  the  other  layer.  The  coherence  length  is  inversely  proportional  to  the  Fermi  wavenumber 
(kp)  and  is  approximately  given  by7 


(4) 

The  critical  current  is  defined  as  the  maximum  interlayer  current  attained  before  self-consistency  is  lost.  After  the 
system  passes  the  interlayer  voltage  corresponding  to  the  critical  current,  the  interlayer  current  drops  sharply  as  the 
system  can  no  longer  sustain  the  condensate  phase  and  interlayer  tunneling  becomes  dominated  by  single  particle 
tunneling7  7 .  We  plot  the  interlayer  current-voltage  characteristic  calculated  for  two  values  of  gate  voltage  in  Fig.  2. 
We  first  examine  the  critical  current  with  a  gate  voltage  of  Vtg  =  —Vbg  (hereafter  referred  to  as  Vc)  =  0.40  V,  a  value 
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FIG.  3.  (color  online)  (a)  Superfluid  fraction  for  spinless  (black  circles,  Nf  =  2)  and  spin-|  (blue  triangles,  Nf  =  4)  symmetric 
double-layer  electron-hole  condensate,  calculated  from  PIMC  simulations  of  40  electron-hole  pairs.  Estimates  of  Tkt  are  made 
where  the  superfluid  fraction  drops  to  1/2,  and  lines  are  a  guide  to  the  eye.  Normalized  historgrams  of  topological  winding 
distributions,  which  measure  superfluid  density,  Eq.  (5),  for  the  (b)  spinless  and  (c)  spin-|  case. 


chosen  to  place  the  system  well  within  the  predicted  regime  for  room-temperature  superfluidity7 .  We  find  that  the 
critical  current  for  the  ballistic  case  is  Ic  «  ISA/iA,  in  good  agreement  with  the  analytic  value  of  I ^  ~  19.5 iiA  while 
IPh  ^  17.9^4  when  the  phonon  interactions  are  included  representing  a  change  of  only  2.7%.  At  Vq  -  0.20  V,  we 
find  a  critical  current  of  Ic  «  10.0 fiA  for  the  ballistic  case  and  I^h  ^3.4  fiA  when  phonon  interactions  are  included,  a 
66.0%  drop  in  critical  current.  This  work  demonstrates  that  if  the  graphene  bilayer  is  operated  at  high  gate  voltages, 
then  the  effects  of  phonon  scattering  can  be  minimized. 


B.  Overall  Scientific  Accomplishments:  Phase  Transitions  with  Increased  Fermion  Flavor 

Interest  in  dipolar  superfluids  has  received  increased  due  in  large  part  to  the  prediction  of  dipolar  superfluid 
behavior  at  or  above  room  temperature  in  double  layer  graphene7-9.  This  is  uniquely  possible  in  graphene  due  to 
the  symmetric  linear  band-structure  and  ability  to  sustain  large  carrier  concentrations  in  two  closely-spaced  layers. 
Yet,  this  prediction  is  not  without  significant  controversy.  As  superfluidity  is  predicted  to  occur  in  the  double  layer 
graphene  system  outside  of  the  quantum  Hall  regime,  additional  fermion  flavors — beyond  the  top  or  bottom  layer 
freedom — may  participate  in  the  phase  transition.  Theoretical  disagreements  over  the  Kosterlitz-Thouless  transition 
temperature  {Tkt)  arise  from  differing  assumptions  about  the  importance  of  these  extra  flavors  for  screening  in 
dipolar  exciton  condensates.  In  the  works  predicting  a  high  transition  temperature  of  Tkt  ~  0.1  Tp  (where  TF  is 
the  system  Fermi  temperature)  the  fermionic  degrees  of  freedom  in  the  system  were  taken  to  be  strongly  correlated 
in  a  gapped  condensate  phase  and  could  not  screen  enough  to  significantly  lower  Tkt-  Other  works  that  predict 
a  low  transition  temperature  of  TKt  ~  10 ~7  Tp  assume  that  screening  from  additional  degrees  of  freedom  add 
independenty  to  effectively  screen  out  the  the  interlayer  interaction10,11.  While  experiment  will  be  the  ultimate 
arbiter  of  the  value  of  Tkt,  many-body  theoretical  approaches  will  play  a  significant  role  in  understanding  the 
nature  of  the  of  the  phase  transition. 

In  Fig.  3  (a),  we  show  our  calculated  superfluid  fraction  as  a  function  of  the  system  temperature.  Our  calculations 
reveal  a  suppressed  Tkt  when  the  additional  spin  degree  of  freedom  is  included.  Superfluid  fraction  is  an  excellent 
way  to  understand  phase  transitions  in  systems  as  it  saturates  at  1  for  low  temperatures  and  asymptotes  to  0  beyond 
the  phase  boundary.  We  denote  the  Tkt  as  the  temperature  where  the  superfluid  fraction  has  dropped  to  0.5  to  be 
consistent  with  previous  work12.  We  see  a  clear  drop  in  TKt  from  TKt  ~  580  K  in  the  spinless  case  to  TKt  ~  300 
K  in  the  spin—  case,  illustrating  the  effect  of  increased  fermion  flavor  on  the  phase  transition.  In  fact,  this  drop  in 
transition  temperature  may  be  understood  without  resorting  to  screening  arguments.  One  would  expect  a  decrease 
in  Tkt  simply  from  reduced  quantum  degeneracy  as  more  fermion  flavors  are  added.  In  an  ideal  Bose  gas,  the 


4 


magnitude  of  the  transition  temperature  is  determined  by  the  condition  n\(TKT )2  ~  1,  where  A (T)  ~  T-1/2  is  the 
thermal  deBroglie  wavelength  and  n  is  the  density  of  identical  Bosons.  As  fermions  pair  into  bosons,  there  must 
be  at  least  Nf  distinguishable  species  of  bosons,  decrease  the  density  n  and  hence,  the  transition  temperature  Tkt 
by  l/Nf.  This  effect  of  reduced  quantum  degeneracy  degeneracy  can  also  be  seen  by  comparing  Tkt  to  the  Fermi 
temperature  (at  our  density,  T^s  =  3086  K  and  Ts/  =  1543  K);  in  each  case  our  estimated  Tkt  is  around  four  or  five 
times  smaller  than  Tp. 

The  calculated  superfluid  fraction  is  estimated  by  the  presence  of  permuting  paths  that  wind  around  the  periodic 
box13.  In  two  dimensions,  the  superfluid  fraction  is  given  by. 
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In  Eq.  (5),  Wx  and  Wy  are  the  topological  windings  (which  are  integer  multiples  of  the  supercell  dimensions)  cor¬ 
responding  to  the  different  path  configurations.  Note  that  for  excitonic  condensates,  we  use  the  number-coupled 
winding,  where  Wx  and  Wy  each  include  a  sum  over  individual  windings  of  electron  and  hole  quasi-particles. 

In  Fig.  5,  the  bottom  two  rows,  (a)  and  (b)  compare  the  exciton  winding  probabilities  for  the  spinless  (a)  and 
spin- 7^  (b)  for  several  temperatures.  At  T  =  150  K,  we  see  similar  behavior  for  both  systems,  as  each  system  is  nearly 
all  superfluid.  Careful  examination  shows  the  histogram  for  the  spinless  case  have  extra  even-odd  structure.  In  an 
exciton,  electrons  pair  with  the  holes  and  wind  as  a  composite  particle  giving  rise  to  peaks  in  the  winding  histogram 
at  even  integers.  Thus  showing  clear  evidence  the  system  is  an  exciton  condensate.  The  even-odd  effect  is  more 
washed  out  in  the  spin—  state,  which  we  attribute  to  the  relatively  small  number,  ten,  of  identical  fermions.  Larger 
simulations  should  sharpen  this  even-odd  signature  of  an  excitonic  condensate.  Note  that  biexcitons  would  have 
strong  peaks  at  multiples  of  four;  we  do  not  see  evidence  of  a  biexcitonic  condensate  in  our  simulations. 

Next  we  examine  real-space  static  correlation  among  the  quasiparticles,  due  to  Coulomb  interactions  and  Pauli 
exclusion.  We  have  collected  the  pair-correlation  functions  g(r),  which  are  easily  calculated  in  PIMC  by  binning  the 
equal-time  pair  distributions  in  a  histogram.  In  contrast  to  other  QMC  methods,  such  as  variational  or  diffusion 
Monte  Carlo,  there  is  no  bias  from  a  trail  wavefunction,  thus  PIMC  gives  essentially  exact  pair  correlations,  aside 
from  errors  associated  with  the  fixed-node  approximation.  In  Fig.  4  we  show  our  calculated  pair  correlation  func¬ 
tions  at  T  =  150  K.  We  find  the  presence  of  an  exchange  hole  in  between  identical  fermions,  which  is  larger  in  the 
Nf  =  4  case  simply  because  the  density  of  particles  of  a  particular  flavor  decreases  as  l/Nf.  We  find  that  there  is  a 
strong,  attractive,  interlayer  correlation  geh,  with  a  peak  at  small  separations  that  is  coincident  with  the  formation  of 
indirectly  bounds  excitons  with  holes  in  the  top  layer  sitting  directly  above  the  electrons  in  the  bottom  layer.  This 
pairing  of  identical  spin  states  in  the  Nf  =  4  case  may  be  an  artifact  of  our  nodal  model,  where  we  have  built  or 
nodal  wavefunction  from  a  Slater  determinant  of  orbitals  pairing  like  spins  only. 

Interestingly,  we  find  a  small  correlation  hole  for  the  Nf  =  4  case,  arising  almost  equally  from  the  electron  and  hole 
of  opposite  spin:  the  particle  in  the  same  layer  (with  the  same  charge)  is  slightly  repelled,  while  the  particle  in  the 
opposite  layer  (with  the  opposite  charge)  is  similarly  attracted.  This  is  a  clear  demonstration  extra  fermion  flavors 
correlate  among  each  other,  even  when  they  do  not  bind  as  condensed  excitons.  This  is  more  dramatically  illustrated 
in  Fig.  4  (c)  and  (d).  In  Fig.  4,  we  plot  the  intralayer  pair  correlation  in  2D  corresponding  to  g and,  as  expected, 
we  see  a  large  correlation  hole  where  Pauli  exclusion  and  Coulomb  repulsion  forces  other  identical  particles  apart. 
However,  when  we  include  spin—  particles  and  examine  the  intralayer  correlation  of  g^  in  Fig.  4  (c),  we  see  a  much 
weaker  correlation  hole.  Note  that  we  again  do  not  see  evidence  of  biexciton  formation,  which  would  appear  as  a 
positive  correlation  of  like-charged,  opposite-spin  particles  in  the  same  layer7 . 

We  do  not  see  a  large  drop  in  Tkt  as  predicted  in  Refs.  10  and  11.  In  an  effort  to  understand  this,  we  explore  the 
dynamic  correlation  functions,  namely  the  polarizability.  In  Refs.  10  and  11,  the  large- Nf  approximation  strongly 
resembles  the  random  phase  approximation  (RPA),  and  hence  the  screening  they  find  is  essentially  the  Lindhard 
function,  with  suppression  of  screening  at  very  small  q  due  a  BCS-like  excitonic  condensate.  In  Fig.  5,  we  show  that 
our  calculated  polarizabilities  are  considerably  different.  To  collect  the  dynamic  correlation  functions  we  sample  the 
polarization  operator  between  different  fermion  flavors  7  and  7', 
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where  nq>7  =  ezq  ri  is  the  density  operator  for  fermion  flavor  7. 

In  Fig.  5  we  show  the  II7  y  (q,  0)  between  different  fermion  flavors  normalized  by  the  2D  density  of  states  as  a 
function  of  wavelength  normalized  by  the  Fermi  wave  vector  at  a  system  temperature  of  T  =  150  K.  Here  we  find 
the  polarization  response  from  identical  particles  in  Nf  =  2  and  Nf  =  4  track  the  Lindhard  function  very  well  in  the 
short  wavelength  limit,  past  q>Tkf,  as  one  would  expect  from  Fermi  liquid  theory. 
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FIG.  4.  (color  online)  Radial  pair  correlation  functions  for  (a)  spinless,  Nf  =  2,  and  (b)  spin-4,  Nf  =  4  cases,  (a)  Intralayer 
exchange-correlation  hole  for  the  same  spin  species,  (d)  Intralayer  correlation  hole  for  opposite  spin  species,  9U(  rw). 


Interestingly,  we  find  similar  behavior  for  both  Nf  =  2  and  Nf  =  4  in  response  to  long  wavelength  excitations 
(q  0).  This  is  expected  as  we  are  clearly  in  the  condensate  regime  in  which  all  carriers  should  be  paired  into 
excitons  and,  therefore,  unable  to  respond  to  perturbations.  The  response  between  particles  in  different  layers  or 
of  different  spin  is  negative,  indicating  that  these  correlations  suppress  screening  from  the  independent  particle, 
RPA  response.  The  response  of  the  exciton  is  quite  strong,  which  the  opposite-spin  flavors  have  a  weak  negative 
response.  Fig.  5(b).  We  have  calculated  the  polarizabilities  at  temperatures  both  above  and  below  T KT,  and  found 
very  little  difference  between  the  different  temperature  traces  for  both  the  Nf  =  2  and  Nf  =  4  cases,  indicating  that 
the  screening  behavior  is  dominated  by  pre-formed  excitons,  with  exist  above  Tkt- 

However,  Fig.  5(c)  shows  that  the  sum  of  inter-  and  intra-layer  polarizabilities  for  Nf  =  4  is  smaller  than  that  of 
Nf  =  2.  While  the  large-N  calculation  finds  enhanced  screening  with  more  fermion  flavors  Nf  which  grows  as  N, 
but  we  find  a  weaker  trend.  The  reason  for  this  difference  is  excitonic  binding  and  correlation  between  different.  In 
the  large-N  expansion,  the  starting  point  is  RPA  like  screening.  Because  each  flavor  can  screen  independently,  it  has 
been  argue  that  the  Coulomb  interaction  because  too  weak  to  form  pairs,  except  for  weak  BCS  pairing  of  states  near 
the  Fermi  surface  when  the  temperature  is  below  Tkt •  In  our  simulations,  our  model  has  strong  excitonic  pairing. 
The  pairing  of  charge  into  neutral  excitons  suppresses  their  ability  to  screen.  This  excitonic  pairing  and  correlation 
between  quasiparticles  of  different  spin  is  not  present  in  the  large- Nf  approximation's  perturbative  expansion. 

This  work  shows  that  the  phase  transition  in  very  closely  separated  layers  of  graphene  is  not  possible  at  room 
temperature.  This  is  because  the  exciton  condensate  is  not  an  effective  screener  of  additional  fermion  flavors  present 
in  the  graphene  system.  That  said,  the  most  interesting  scientific  breakthrough  this  year  has  been  the  successful 
implementation  of  linear  response  theory  within  PIMC.  It  has  given  us  a  tool  to  examine  very  interesting  and  ex¬ 
perimentally  relevant  phenomena  without  needing  a  large  number  of  approximations  and  with  full  many-body 
interactions. 


C.  Overall  Scientific  Accomplishments:  Transport  in  Clean  and  Disordered  Zero-Field  Semiconductor  Bilayer  Systems 

For  almost  two  decades  coupled  semiconductor  quantum  well  systems  have  revealed  interesting  effects  deep 
within  the  quantum  Hall  regime.  Comparatively  speaking,  there  have  been  relatively  few  studies  of  the  bilayer  sys¬ 
tem  outside  of  the  quantum  Hall  regime.  Here  we  do  not  expect  to  see  a  phase  transition,  but  the  system  may  still 
exhibit  transport  properties  similar  to  those  characteristic  of  the  quantum  Hall  pseudospin  ferromagnets,  such  as 
an  enhancement  on  interlayer  tunneling  conductance14.  Previous  mean-field  calculations  show  that  semiconductor 
bilayers  may  exhibit  interlayer  exchange  enhancement  characterized  by  enhancement  factor  S,  but  the  value  of  S 
is  moderate15.  Following  their  work,  we  examine  the  interlayer  transport  properties  of  both  clean  and  disordered 
semiconductor  bilayer  systems  without  magnetic  fields.  We  find  that  transport  characteristics  similar  to  quantum 
Hall  ferromagnets  may  persist  outside  of  the  quantum  Hall  regime  and  we  discuss  how  the  the  physical  properties 
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FIG.  5.  (a)  Interlayer  and  intralayer  polarizabilities  for  Nf  —  2  versus  the  perturbation  wavevector  normalized  by  the  Fermi  wave 
vector,  (b)  Interlayer  and  intralayer  polarizabilities  for  Nf  =  4  normalized  by  the  Fermi  wave  vector,  (c)  Sum  of  the  interlayer 
and  intralayer  polarizabilites  for  Nf  —  2  and  Nf  =  4  normalized  by  the  Fermi  wave  vector.  In  each  plot  the  polarizabilities  are 
plotted  against  the  Lindhard  response  function  (grey  line)  in  2D  and  have  been  normalized  by  the  density  of  states  at  the  Fermi 
energy  in  2D,  No . 
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FIG.  6.  Device  cross-section  in  x,z  direction.  Top  and  bottom  Alo.9Gao.1As  barrier  thickness  is  60  nm,  GaAs  quantum  well 
thickness  is  15  nm,  and  Alo.9Gao.1As  barrier  thickness  is  1  nm.  A  channel  is  1.2  iim  long  and  7.5  (im  wide.  An  assumed  electron 
density  in  2  dimensional  electron  gas  (2DEG)  is  2.0  x  10locra-2. 


of  these  exchange  enhanced  transport  signatures  change  as  we  include  material  disorder.  This  work  is  of  particular 
importance  as  it  provides  us  with  a  system  which  is  both  experimentally  realizable  and  readily  accessible  to  exper¬ 
imentalists.  We  study  this  system  to  provide  a  pathway  to  understand  the  validity  of  our  graphene  theory  in  the 
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FIG.  7.  Plot  of  the  differential  conductance  of  the  bilayer  system  at  S  =  4  (orange  triangle),  S  =  2  (blue  circle),  and  S  =  1  (opened 
rectangular)  as  a  function  of  interlayer  bias  Vb  —  Vt  at  OK.  In  the  inset,  the  same  plot  in  negative  interlayer  bias  with  same  y 
axis  is  plotted  with  x  axis  in  log  scale.  In  small  bias  window,  almost  constant  interlayer  transconductance  is  observed.  After  the 
interlayer  bias  reaches  Vint  ~  10 fiV,  an  abrupt  drop  in  transconductance  can  be  seen  in  the  interacting  electron  cases  of  S  =  2, 4. 
We  plot  the  differential  conductance  on  a  logarithmic  scale  in  the  inset  figure  to  clearly  illustrate  the  observed  drop  in  interlayer 
conductivity. 


absence  of  graphene  experiments  and  with  the  effects  of  screening  making  the  observation  of  the  condensate  phase 
doubtful. 

A  schematic  feature  of  our  device  is  presented  in  Fig.  6.  In  this  work,  we  consider  for  an  Alo.9Gao.1As/GaAs 
bilayer  heterostructure  with  60  nm  top  and  bottom  Alo.9Gao.1As  barriers  which  act  to  isolated  the  coupled  quan¬ 
tum  wells  from  electrostatic  gates.  Our  system  is  comprised  of  two  15nm  deep  GaAs  quantum  wells  with  a  1  nm 
Alo.9Gao.1As  barrier  between  the  two  quantum  wells.  Since  bilayer  experiments  have  taken  place  at  very  low  tem¬ 
perature  around  a  few  tens  of  mK,  we  may  write  the  Hamiltonian  in  the  form  of  decoupled  ID  longitudinal  channels 
in  transport  (x)  direction  taking  account  of  eigenenergy  of  transverse  ( y )  direction  properly16.  The  system  is  coupled 
with  a  3D  Poisson  solver  to  correctly  account  for  the  electrostatics.  Bottom  contacts  are  biased  and  current  is  ex¬ 
tracted  out  from  top  contacts  ( Vbl  =  Vbr  =  Vint ,  Vtl  =  Vtr  =  0).  From  our  analysis,  the  critical  current  of 
this  system  is  determined  by  single  particle  tunneling  amplitude  and  enhancement  factor  S.  The  resulting  critical 
currents  are  Ic  =  1.80,  4.52  nA  for  S  =  2,4  respectively  which  show  a  good  agreement  with  mean-field  calculation 
results,  =  2.14,  4.00  nA. 

To  get  a  clear  picture  of  how  interlayer  interactions  modify  the  resultant  interlayer  transfer  characteristics,  we 
force  the  system  out  of  the  linear  response  regime.  In  Fig.  (7)  ,  we  compare  the  interlayer  transport  result  with 
the  non-interacting  case  (S  =  1).  We  see  that  all  of  the  curves  are  sharply  peaked  near  zero  bias,  more  strongly  in 
the  case  of  the  inter  acting-system  simulations  (S  =  2,4)  which  is  expected.  After  the  critical  current,  we  see  that 
the  interlayer  conductivity  drops  appreciably,  as  is  the  case  in  graphene  or  the  quantum  Hall  system.  However,  the 
important  thing  to  point  out  here  is  that  the  semiconductor  bilayer  system  has  a  critical  current  which  does  not  occur 
when  the  in-plane  pseudospin  phase  angle  reaches  | .  This  is  due  to  the  charge  imbalance  which  develops  between 
the  layers  causing  the  out-of-plane  pseudospin  field  to  precess  faster  separating  the  layers  energetically  with  smaller 
applied  interlayer  bias. 

With  our  understanding  of  the  system  in  clean  limit,  we  are  interested  in  AlGaAs/GaAs  system  with  impurities 
which  is  naturally  occurred  during  a  fabrication.  In  order  to  simulate  device  characteristics,  we  adopt  a  model  of 
ionized  impurities  and  roughness  at  the  interface  of  heterojunction  in  AlGaAs/GaAs  system17.  Random  disorders 
can  be  inserted  into  the  system  Hamiltonian  by  adding  uniform  distribution  having  a  window  defined  by  mean  free 
path,  A.  The  simulation  is  carried  out  with  NEGF  formalism  in  two  dimensions  coupled  with  3D  Poisson  solver.  A 
mean  free  path  we  examine  in  this  study  adds  a  significant  amount  of  perturbation  into  onsite  energy  up  to  ~  20%  of 
the  self-consistent  Hartree  potential  at  clean  limit.  The  perturbation  in  turn  affects  the  expectation  values  of  the  in¬ 
plane  and  out  of  plane  pseudospin  fields  and  alters  transport  characteristics  of  the  system.  Clearly,  disorder  causes 
some  localization  in  the  transport  (x)  directional  pseudospin  fields  and  pseudospin  angles  of  transport  orbital  are 
no  longer  well  aligned  and  locally  precessed  even  at  the  low  bias  window.  The  overall  effect  of  disorder  results  in 
enhancement  on  interlayer  tunneling,  but  less  than  the  enhancement  in  clean  limit.  The  random  scatterers  are  added 


8 


FIG.  8.  Plot  of  differential  conductance  as  a  function  of  interlayer  bias.  We  plot  clean  limit  of  assumed  enhancement  factor  S  =  2 
with  non-interaction  case  S  m  l.  We  also  present  a  result  with  a  disorder  level  of  A  =  0.3  gm  with  error  bars  from  four  different 
disorder  configurations.  Simulation  conditions:  Interlayer  bias  is  only  applied  at  Vbr  =  Vint ,  and  rest  of  the  contacts  is  zero. 
An  assumed  2DEG  density  is  ne  =  3.0  x  10“ 10 cm~ 2  and  quantum  well  has  Length  =  Width  =  1.2  fim. 


up  to  A  =  0.3  iim  and  the  end  result  is  a  reduction  of  the  interlayer  conductance  to  its  non-interacting  value  as  shown 
in  figure  8.  The  second  major  effect  of  the  disorder  is  broadening  of  tunneling  peak  described  by  T  =  h/r  where  r 
is  quantum  life  time18.  The  predicted  broadening  (T/2)  for  A  =  0.3  iim  is  around  82  /neV  and  corresponding  result 
from  NEGF  calculation  lies  roughly  around  80  /ieV  in  Fig.  8.  This  indicates  that  the  interlayer  transport  is  no  longer 
controlled  by  single  particle  tunneling  amplitude  and  instead  is  controlled  by  the  material  disorder. 

This  work  gives  an  insight  into  electron-electron  bilayer  system  under  zero-magnetic  field,  a  system  which  is 
experimentally  realizable..  We  find  that  interlayer  transport  in  this  system  is  limited  by  the  strength  of  correlation 
between  two  layers  and  charge  imbalance  of  the  system.  We  find  localized  pseudospin  field  with  disorder  and  their 
net  effect  suppresses  coherent  tunneling  transport.  With  significant  amount  of  disorder,  interlayer  transport  tends  to 
be  dependent  on  material  disorder. 


D.  Final  Scientific  Accomplishments:  Dipolar  Superfluids  in  Disordered  Systems 

Our  simulations  in  the  ideal  system  show  two  different  regions  in  the  graphene  bilayer  -  a  region  near  the  contacts 
where  injected  carriers  trigger  a  crossed  Andreev  Reflection  event  and  the  channel  region  where  all  transport  is 
attributed  to  superfluid  current.  When  an  electron  is  injected  with  energy  lower  than  the  superfluid  gap,  it  will 
trigger  an  Andreev  Reflection  similar  to  what  happens  in  a  superconductor.  In  the  exciton  superfluid  case,  however, 
an  injected  electron  will  drive  an  exciton  across  the  channel,  with  a  retroreflected  hole  coming  out  the  bottom  left 
contact  to  conserve  momentum.  The  effect  of  the  crossed  Andreev  reflection  is  thus  a  superfluid  intralayer  current 
and  the  quasiparticle  interlayer  current.  The  quasiparticle  current  evanescently  decays  into  the  channel,  penetrating 
a  certain  length  into  the  superfluid  region.  This  distance,  known  as  the  coherence  length,  is  found  to  be  Lc  =  5  nm. 
We  use  this  observation  to  break  down  the  disordered  case  study  into  two  scenarios  -  where  disorder  is  located  only 
in  the  coherence  length  (Contact  Disorder)  and  where  disorder  is  located  away  from  the  coherence  length  (Channel 
Disorder). 

The  main  results  of  our  work  are  shown  in  Fig.  9.  At  a  vacancy  concentration  of  0%,  the  critical  current  is  18.4  fiA. 
To  the  left  of  the  ideal  point  we  plot  the  effect  of  channel  disorder  on  critical  tunneling  current.  We  find  that  vacancies 
induce  a  local  density  of  states  (LDOS)  in  graphene  layers  that  bring  about  occupied  states  within  the  superfluid  gap 
and  act  to  locally  lower  the  superfluid  density,  mexc.  For  small  vacancy  concentrations  (1%  vacancies  placed  ran¬ 
domly  throughout  the  channel),  this  has  little  effect  on  mexc  because  the  disorder  is  sufficiently  isolated.  However, 
as  the  channel  is  made  more  disordered,  superfluid  density  begins  to  drop  non  linearly  in  the  channel.  The  average 
superfluid  density  saturates,  however,  due  to  the  unchanging  contribution  from  the  ideal  contact  regions,  decreasing 
to  roughly  70%  of  its  value  at  a  4%  vacancy  concentration.  Because  critical  tunneling  current  is  proportional  to  mexc , 
we  find  good  agreement  with  this  explanation  via  the  approximate  square-root  dependence  seen  in  Fig.  9.  The  error 
bars  show  that  the  configuration  of  the  disorder  has  little  effect  on  Ic/  as  average  superfluid  density  remains  robust 
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FIG.  9.  Statistical  variation  of  critical  interlayer  current  for  several  runs  at  the  bias  just  before  the  condensate  is  broken.  Vacancies 
which  are  more  than  a  distance  of  one  coherence  length  away  from  the  top  layer  contacts  are  to  the  left  of  the  ideal  (0%)  case, 
depicted  by  a  black  circle,  while  vacancies  which  are  within  one  coherence  length  from  the  top  layer  contacts  are  to  the  right.  The 
excitonic  condensate  is  lost  in  both  cases  for  top  layer  vacancy  concentrations  larger  than  4%.  The  solid  lines  represent  analytic 
calculations  for  Ic  with  modified  values  of  superfluid  density  (blue  line)  and  width  (red  line). 


regardless  of  the  disorder's  random  configuration.  In  both  cases,  no  self-consistent  solution  was  found  for  vacancy 
concentrations  greater  than  4%. 

We  see  a  clear  linear  dependence  in  the  contact  disordered  regime.  Injected  carriers  become  increasingly  likely  to 
be  backscattered  by  the  impurities  near  the  contacts  before  Andreev  reflection  can  occur.  This  results  in  a  significant 
reduction  in  the  area  in  which  quasiparticle  tunneling  can  occur  near  the  contact,  reducing  interlayer  conductivity 
significantly.  As  critical  tunneling  current  is  linearly  dependent  on  the  width  of  the  system,  this  provides  a  simple 
explanation  for  the  physics  of  the  linear  decrease  in  Ic.  Assuming  that  each  vacancy  removes  0.5  nm  of  effective 
width  from  the  system,  we  plot  the  resulting  relationship  in  Fig.  9  with  a  red  line  denoting  I ^  with  a  decreasing 
width.  Deviation  in  Ic/  as  seen  by  the  large  error  bars  in  the  contact  disordered  case,  is  attributed  to  the  non  linear 
randomized  bunching  of  disorder. 

This  work19  provides  a  comprehensive  study  on  the  effect  of  short-range  scatterers  in  double  layer  graphene  sys¬ 
tems  that  must  be  taken  into  account  for  this  system  to  be  experimentally  realized.  Unfortunately,  room  temperature 
condensation  is  subject  not  only  to  defects  in  fabrication,  but  also  carrier-phonon  interactions. 


E.  Overall  Scientific  Accomplishments:  Dipolar  Superfluids  in  3D  Time-Reversal  Invariant  Topological  Insulators 

Spatially  segregated  monolayers  of  graphene  have  been  both  theoretically7-9  and  experimentally20  explored  for 
signatures  of  excitonic  superfluidity.  While  signs  of  interlayer  correlation  are  experimentally  observed,  additional 
fermionic  degrees  of  freedom,  or  flavors,  screen  the  strength  of  the  interlayer  interaction21  making  the  observation 
of  dipolar  exciton  superfluid  in  graphene  multilayers  challenging. 

The  advent  of  time-reversal  invariant  topological  insulators  (TI)22,23  has  brought  renewed  interest  in  finding  DES 
in  condensed  matter  systems.  In  sufficiently  thin  TI  films,  it  has  been  proposed  that  spatially  segregated  surface 
electrons  and  holes  may  bind  into  a  topological  dipolar  intersurface  exciton  superfluid  (TDIES).  To  this  point,  exist¬ 
ing  approaches  to  TDIES  have  considered  strictly  two-dimensional  Dirac  surface  states  separated  by  an  insulating 
spacer24-28.  Yet  the  existence  of  a  TDIES  in  three-dimensions  is  not  a  foregone  conclusion  based  on  two-dimensional 
surface  state  analysis.  The  most  obvious  drawback  being  that  in  a  3D  TI,  each  of  the  surfaces  is  interconnected  and 
there  exists  no  obvious  mechanism  to  segregate  the  electron  and  hole  layers,  as  in  other  proposed  systems. 

We  begin  in  Fig.  10,  where  we  schematically  show  the  system  we  consider.  We  apply  top  and  bottom  potentials  of 
opposite  polarity  to  induce  electrons  on  the  top  surface  and  holes  on  the  bottom  surface.  We  attach  contacts  Cl  -  C4  to 
the  top  and  bottom  surfaces  on  the  left  and  right  sides  of  the  TI  through  which  current  may  be  injected  and  extracted 
as  seen  in  Fig.  10.  As  for  the  major  findings  resulting  from  this  work,  we  find  that  the  proper  exciton  pairing  order 
parameter  is  purely  imaginary  which  is  necessary  to  properly  account  for  the  system  dynamics.  Furthermore,  we 
find  that  the  exciton  order  parameter  is  p-wave  and  that  electrons  will  only  bind  with  holes  of  the  same  chirality 
from  different  orbitals.  We  find  that  TDIES  in  time-reversal  invariant  TI  thin  films  prefer  to  bind  into  excitons  with  a 
dominant  many-body  energy  which  prevents  single  particle  intersurface  transport  as  the  intersurface  polarizablity 
drops  to  zero.  This  allows  for  the  observation  of  superfluid  behavior  out  of  the  quantum  Hall  regime  and  without 
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FIG.  10.  Schematic  of  topological  insulator  thin-film  system  under  consideration.  The  top  and  bottom  surfaces  are  assumed  to 
contain  equal  numbers  of  electron  and  holes,  respectively. 


FIG.  11.  A  transmission  between  Cl,  C2,  C3,  and  C4  of  the  system  (a)  without  and  (b)  with  intersurface  interaction  at  a  bias 
of  Vbias  —  5  mV.  (a):  When  there  is  no  interaction  and  the  system  is  not  gapped,  current  flows  across  the  device  resulting  in 
T12  (or  Z34)  dominating  the  transport,  (b):  When  the  system  is  gapped,  however,  the  quasiparticle  undergoes  a  similar  process 
to  Andreev  reflections,  resulting  in  large  T13  and  T24  as  current  flows  from  top  surface  to  bottom  surface  with  no  transmission 
across  either  of  the  top  Ti2  or  bottom  surfaces  Z34. 


the  need  to  artificially  segregate  the  surfaces.  Finally,  we  find  that  the  presence  of  a  topological  superfluid  may  be 
electrically  detected  by  the  presence  of  equal  and  opposite  intersurface  contact  currents,  as  can  best  be  seen  from  the 
individual  transmission  in  Fig.  11. 


III.  FINAL  CONCLUSIONS  AND  FUTURE  OUTLOOK 


The  goals  of  this  3  year  project  have  been  to  (i)  analyze  the  physics  of  high-temperature  superfluidity  in  bilayer 
graphene  systems  with  a  goal  of  understanding  its  role  in  post-CMOS  architecture  and  (ii)  develop  the  path  integral 
quantum  Monte  Carlo  technique  as  a  new  flexible  approach  for  many-body  quantum  systems  at  finite  temperature. 
With  regards  to  the  first  goal  of  determining  the  role  of  segregated  layers  of  tightly  spaced  graphene  for  post-CMOS 
architectures,  the  conclusions  of  our  work  show  that  the  bilayer  graphene  system  has  transport  properties  in  the 
superfluid  phase  which  outperform  other  proposed  post-CMOS  devices.  Yet  our  path  integral  quantum  Monte  Carlo 
calculations  show  that  the  superfluid  phase  transition  is  occluded  by  additional  fermionic  degrees  of  freedom  and 
this  effectively  lowers  the  transition  temperature  to  well  below  room  temperature.  We  have  shown  that  topological 
insulators  provide  a  much  more  promising  system  for  high-temperature  device  operation  and  these  systems  should 
be  further  explored.  As  to  the  second  goal  of  this  project,  the  development  of  efficient  path-integral  quantum  Monte 
Carlo  algorithms  for  finite  temperature  many-body  simulations,  we  have  successfully  pushed  this  technique  farther 
than  ever  before.  In  the  sections  outlined  above,  we  have  shown  the  ability  to  collect  static  and  dynamic  response 
functions  for  fermionic  many-body  systems  to  elucidate  the  superfluid  transition  temperature.  This  represents  a 
powerful  new  technique  from  which  the  full  interacting  conductivity  may  be  extracted.  I  feel  that  the  best  way 
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forward  with  regards  to  the  path-integrals  would  be  to  take  a  step  back  from  applying  them  to  emergent  systems, 
whose  responses  are  not  known,  and  calibrate  the  technique  against  other  well-known  systems  such  as  Luttinger 
liquids  in  quantum  wires. 


1  J.  Cayssol,  B.  Huard,  and  D.  Goldhaber-Gordon,  Phys.  Rev.  B  79,  075428  (2009). 

2  J.  Su  and  A.  H.  MacDonald,  Nat.  Phys.  4,  799  (2008). 

3  S.  Datta,  Quantum  transport :  atom  to  transistor  (Cambridge  University  Press,  Cambridge  UK  ;;New  York,  2005). 

4  E.  H.  Hwang  and  S.  Das  Sarma,  Physical  Review  B  77, 115449  (2008). 

5  K.  M.  Borysenko,  J.  T.  Mullen,  E.  A.  Barry,  S.  Paul,  Y.  G.  Semenov,  J.  M.  Zavada,  M.  B.  Nardelli,  and  K.  W.  Kim,  Phys.  Rev.  B  81, 
121412  (2010). 

6  Y.  Yoon,  D.  E.  Nikonov,  and  S.  Salahuddin,  Applied  Physics  Letters  98,  203503  (2011). 

7  H.  K.  Min,  R.  Bistritzer,  J.  J.  Su,  and  A.  H.  MacDonald,  Phys.  Rev.  B  78, 121401  (2008). 

8  M.  J.  Gilbert  and  J.  Shumway,  J.  Comput.  Electron.  8,  51  (2009). 

9  C.  H.  Zhang  and  Y.  N.  Joglekar,  Phys.  Rev.  B  77,  233405  (2008). 

10  M.  Y.  Kharitonov  and  K.  B.  Efetov,  Phys.  Rev.  B  78,  241401R  (2008). 

11  M.  Y.  Kharitonov  and  K.  B.  Efetov,  Semicond.  Sci.  Technol.  25,  034004  (2010). 

12  E.  L.  Pollock  and  D.  M.  Ceperley,  Phys.  Rev.  B  36,  8343  (1987). 

13  D.  M.  Ceperley,  Rev.  Mod.  Phys.  67,  279  (1995). 

14  I.  B.  Spielman,  J.  P.  Eisenstein,  L.  N.  Pfeiffer,  and  K.  W.  West,  Phys.  Rev.  Lett.  84,  5808  (2000). 

15  L.  S-acutewierkowski  and  A.  H.  MacDonald,  Phys.  Rev.  B  55,  R16017  (1997). 

16  R.  Venugopal,  Z.  Ren,  S.  Datta,  M.  S.  Lundstrom,  and  D.  Jovanovic,  Journal  of  Applied  Physics  92,  3730  (2002). 

17  T.  Ando  and  H.  Tamura,  Phys.  Rev.  B  46,  2332  (1992). 

18  N.  Turner,  J.  T.  Nicholls,  E.  H.  Linfield,  K.  M.  Brown,  G.  A.  C.  Jones,  and  D.  A.  Ritchie,  Phys.  Rev.  B  54, 10614  (1996). 

19  B.  Dellabetta  and  M.  J.  Gilbert,  Journal  of  Physics:  Condensed  Matter  23,  345302  (2011). 

20  S.  Kim,  I.  Jo,  J.  Nah,  Z.  Yao,  S.  K.  Banerjee,  and  E.  Tutuc,  Phys.  Rev.  B  83, 161401  (2011). 

21  J.  Shumway  and  M.  J.  Gilbert,  Phys.  Rev.  B  85,  033103  (2012). 

22  M.  Z.  Hasan  and  C.  L.  Kane,  Rev.  Mod.  Phys.  82,  3045  (2010). 

23  X.  L.  Qi  and  S.-C.  Zhang,  Rev.  Mod.  Phys.  83, 1057  (2011). 

24  B.  Seradjeh,  J.  E.  Moore,  and  M.  Franz,  Phys.  Rev.  Lett.  103,  066402  (2009). 

25  G.  Y.  Cho  and  J.  E.  Moore,  Phys.  Rev.  B  84,  165101  (2011). 

26  D.  Tilahun,  B.  Lee,  E.  M.  Hankiewicz,  and  A.  H.  MacDonald,  Phys.  Rev.  Lett.  107,  246401  (2011). 

27  Z.  Wang,  N.  Hao,  Z.-G.  Fu,  and  P.  Zhang,  ArXivil  106.5838. 

28  E.  G.  Moon  and  C.  Xu,  EPL  97,  66008  (2012). 


